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' Abstract. In quantum Bayesian inference problems, any conclusions drawn from a finite number 

^ , of measurements depend not only on the outcomes of the measurements but also on a prior. Here 

■ we show that, in general, the prior remains important even in the limit of an infinite number of 
Q>^ I measurements. We illustrate this point with several examples where two priors lead to very different 

conclusions given the same measurement data. 
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§ ■ INTRODUCTION 

^ : 

I— I, This paper addresses the problem of inference in quantum mechanics in a very general 

I I setting. Consider a sequence of quantum systems, each with Hilbert space J^f which, for 

>■ ■ simplicity, is assumed to have finite dimension. Now some of the systems are measured. 

'nI" . What conclusions can be drawn from the measurement outcomes? 

It is often useful to consider less general situations. One of the most studied quantum 
inference problems is quantum state estimation, which is frequently described as follows 
\^ . dll]. Each system is assumed to have the same unknown state a, where cr is a density 

O I operator on Jf. A sequence of appropriately chosen measurements is then used to deter- 

mine o. The most general single-system measurement is described by a POVM, which 
is a collection of positive operators, {^i, . . . ,Er}, acting on J^, such that Y^k^k is the 
. ^ '. identity operator. The index k E { 1 , . . . , r} labels the possible measurement outcomes. If 

^ I the state of the system is p, the probability of obtaining outcome k is given by tr(p£',t)- 

■ The Bayesian solution to the quantum state estimation problem is straightforward. Its 
starting point is a prior probability distribution Pprior(p) on the set of density operators 
on J^. Assume that the POVM {Ej^} is measured for a single system. If the outcome of 
the measurement is k, the distribution Pprior(p) is updated using a quantum version of 
Bayes's rule resulting in the posterior distribution 

This process is iterated, each time using a fresh copy of the system and possibly a dif- 
ferent measurement, and each time setting Pprior(p) equal to the previously obtained 
Pposterior(p)- Given some mild conditions on the original prior distribution, and assum- 
ing that the measurements are appropriately chosen, it can be shown that the iterated dis- 
tribution approaches a delta distribution centered at the unknown state o. For instance. 



0^ 
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this convergence holds if the original Pprior(p) is nonzero for all p and the measured 
POVM is informationally complete [4]. 

In practice, the iteration will have to stop after a finite number of measurements. 
The posterior distribution will then generally depend on the initial prior distribution. 
From a Bayesian perspective, this dependence on the prior distribution is unavoidable. 
Finding the appropriate mathematical form of the prior distribution is a central part of 
the Bayesian approach to quantum state estimation. 

There are a number of (non-Bayesian) state estimation methods that attempt to cir- 
cumvent the dependence of the estimate on the prior distribution [1]. An important ex- 
ample is maximum likelihood estimation. We are not going to discuss non-Bayesian 
methods any further in this paper. 

In the next section, we will show that quantum state estimation as described above is a 
special case of a far more general quantum inference scenario in which priors retain their 
importance even in the limit of an infinite number of measurements. This is followed by 
a section illustrating this point with some examples. The paper concludes with a brief 
discussion. 



GENERAL QUANTUM INFERENCE PROBLEM 

There is something peculiar about the quantum Bayes rule ([T]). Since a density operator 
p encodes the probabilities for outcomes of quantum measurements, the rule mixes 
two kinds of probabilities, namely the "classical" probabilities Pprior and ^posterior on the 
one hand, and the "quantum" probabilities encoded in p on the other hand. We will now 
show that the rule Q is a special case of a more general rule phrased entirely in terms of 
quantum states. It is therefore not necessary to make the distinction between two kinds 
of probability. 

There is a general updating rule built into quantum theory. To describe a quantum 
measurement fully, one needs to provide, in addition to the POVM {^i, . . . ,E,-}, a set of 
Kraus operators Aj^j such that = LjAytjAjy for = 1, . . . ,r [5]. If the measurement 
gives outcome k for a system in state p, the state of the system after the measurement 
will be 

Pk = tv{pEk)-'Y,AkjpAlj. (2) 

j 

This rule has an interpretation very similar to the Bayes rule. If we call p the prior state 
and p^ the posterior state, we see that the prior state is changed into the posterior state 
upon acquisition of the data k. 

To see that the quantum Bayes rule ([1]) is a special case of the Kraus rule Q 10], 
consider a sequence of quantum systems each with Hilbert space as before. We 

(n) 

define a prior on our sequence of systems as a sequence of n-system states p^^^^^, 

in) 

n= 1,2,..., where each Pp^/oj. is a density operator on the n-fold tensor-product Hilbert 

space = ® ■ • • ® jr, and where p^^^^ = ^^n+iPpZ!^ for all n > 1 . By tr„+i we 
denote the trace over the {n + l)-th system. In words, each member of the sequence is 
obtained from the next by tracing over the additional system. 



Now assume that the first system is measured and an outcome k is obtained. Applying 
the rule ^ we find, for any n > 1, that the state of the first n systems after the 
measurement is 

Pl"^=tr(P^l£.)-T^^.Ol' (3) 

j 



where it is understood that the operators Aj^j and E]^ act on the first system only. By 
tracing over the first system, we obtain what we call the posterior on our sequence of 
systems, 

pSLor = triPr^^ (n = l,2,...). (4) 

The posterior is again a sequence of states; its n-ih. member is obtained from the (« + 1)- 
th member of the prior by measuring and then discarding the first system. The posterior 

has the property Pposjerior = t^«+iPposterior all n > 1 and thus has the form of a prior. 
We can therefore iterate the above procedure, each time setting the prior equal to the 
posterior obtained in the previous iteration. 

Given the prior and the measurement data, the posterior is the unique sequence of 
states for the remaining (i.e., not yet measured) systems. In this sense the posterior 
constitutes the unique correct solution of the quantum inference problem. In particular, 
the one-system state Pposterior marginal state for the first unmeasured system. This 
state is sometimes called the Bayesian mean estimator. 

To recover the familiar rule ([I]), only one simple additional assumption has to made, 

in) 

namely that for any n > \ the state p„^-' is symmetric under permutations of the 

J — I 

n systems. In this case we say that the prior is exchangeable [|6[[. Given this extra 
assumption, it is the content of the quantum de Finetti theorem [0, @] that the prior 
can be written as 

Ppnor = / dpp^noAp)p®" (n= 1,2,...) , (5) 

where Pprior(p) is a probability distribution on the space of single-system density oper- 
ators, and p®" is the n-fold tensor product p • ■ ■ C?) p. It is not difficult to establish [0] 
that a measurement on the first system with outcome k will lead to the posterior 

PpoLnor = /^Pi^postenor(p)p®" (n=l,2,...), (6) 

with ;?posterior(p) given by the rule ©• 

As was pointed out in the introduction, in the limit of an infinite number of iterations, 
/'posterior (p) typically approaches a delta function which is independent of the detailed 
functional form of Pprior(p)- This conclusion, however, depends crucially on the assump- 
tion that the prior is of the form ([5]), i.e., that the prior is exchangeable. In other words, 
even in the limit of an infinite number of iterations, conclusions depend on the prior. All 
we can say is that some details of the prior become irrelevant in this limit. 

Exchangeable priors are an important class of priors that are used so frequently that it 
is sometimes overlooked that exchangeability is an assumption. Making this assumption 
is equivalent to choosing a prior from a restricted set. What we have therefore shown in 
this section is that conclusions drawn in Bayesian quantum inference situations generally 



depend on the prior as well as measurement data, even in the limit of an infinite number 
of measurements. In the next section, we illustrate this point with some examples. 



EXAMPLE PRIORS 



The first example is a sequence of qubits, i.e., two-dimensional quantum systems. We 
denote by |0) and 1 1) two orthogonal basis states and consider three priors. 
Our first prior is exchangeable and given by 



pi")= jdw)pa{w)){w){^f\r'' ("=1,2,...) 



(7) 



where d\'^f) Pa{\^f)) is the Haar measure on the space of pure one-qubit states. 

Our second prior is also exchangeable, but consists of a sequence of pure product 
states. It could be called a Rosenkrantz and Guildenstern prior [7]. This is the state 
one might assign to a quantum random number generator manufactured by a trusted 
company. It is given by 



Pb 



(n) 



^(|0) + |1))((0| + (1|))®" (n=l,2,...). 



(8) 



Our third prior is not exchangeable. We call it a counter-inductive prior [8] because 
it leads one to predict outcomes that are the opposite of what an argument by induction 
would suggest. In particular, we will see that updating this prior after a string of m 
identical measurement outcomes, the probability for obtaining the opposite outcome in 
the next measurement approaches 1 as m increases. The counter-inductive prior is given 
by 



£2-^( |0)(0|®^®|i)(i|®("-*) + |i)(i|®^®|0)(0|®("-^) 

k=l 



+ £2-^^(|0)(0|®« + |l)(l|®"jJ (n=l,2,...), 

k=n 

where the normalization constant JV is determined by the equation 

CO 

1 = 2^ £ 2" 

k^\ 



(9) 



(10) 



It is not difficult to check that this sequence of states satisfies the defining condition of a 



(n+l) c ^1 

tr„+iPc ' forn > 1. 



prior, pi 

We now imagine that a sequence of von Neumann measurements in the {|0),|1)} 
basis is carried out and that each measurement produces the same outcome, 0. The 
measurement data thus consist of a string of zeros. These data are used to update 
iteratively each of our three priors above. For each prior, we compute the marginal one- 
system posteriors in the limit of an infinite number of iterations. 



For the two exchangeable priors, we obtain easily the limits 

pa ^ ^\0){0\ (number of iterations ^ oo) (11) 

and 

p^^^ ^ -(|0) + |1))((0| + (1|) (number of iterations ^oo). (12) 
To compute the limit for the counter-inductive prior Pc , we first compute the posterior 

(n) 

after m iterations, Pc,m- We find, for n = 1,2,..., 

, , . W-l oo 

p^;'^ = ^ M 1 ) ( 1 1 8" + £ 2-^(2m+/:) | Q) (0 1 ® 1 1 ) ( 1 1 ® + 2-^(2m+/:) | Q) (0 1 ®" j , 

k=l k=n 

(13) 

where .A^ is determined by the condition 

oo 

k=0 

implying that lim A^ = 1. Hence 

W— »oo 

pi^^ — > 1 1) (1 1 (number of iterations — > oo). (15) 

Clearly, the three priors lead to radically different conclusions for the same infinite 
sequence of data. 

We now move on to our second example, where we compare two priors, again for a 
sequence of qubits. For our first prior, we choose a generic exchangeable prior, 

pf = ldpp,{p)p'^" (n=l,2,...), (16) 

where dp pd{p) is a measure on one-qubit density operators, i.e., density operators on a 

two-dimensional Hilbert space, J^. We assume that Pd{p) is nonzero for all p on M". 
This prior entails that there is no entanglement between the qubits. 

By contrast, our second prior, though identical with our first on the single system 
marginals, does not rule out entanglement between pairs of qubits. For even numbers of 
systems, it is defined by 

pF") = I Jc7p,(a) a®" (n=l,2,...), (17) 

where do pe{<y) is a measure on the space of two-qubit density operators, i.e., density 
operators on the four-dimensional Hilbert space J^f (8) J^. For odd numbers of systems, 
the prior is defined by 

pf = tr2„pf (n=l,2,...). (18) 



Assume now that a sequence of informationally complete measurements is performed 
on the sequence of qubits. As before, the measurement data are used to update both our 
priors iteratively. We assume that the data are such that for the second prior, the marginal 
two-system posterior converges to an entangled two-qubit state, e.g., the maximally 
entangled state puE = ^ ( 1 00) + 1 1 1 ) ) ( (00 1 + ( 1 1 1 ) . We thus assume that 

(2) 

pe Pme (number of iterations ^ oo). (19) 

Since for the maximally entangled state, both marginal states are equal to the totally 
mixed state pM = ^(|0)(0| + |1)(1|), it follows that given the same data, the marginal 
one-system posterior for our first prior converges to pM, 

pI^^ pM (number of iterations oo), (20) 
and the marginal two-system posterior converges to 

(2) 

py ^ Pm^Pm (number of iterations ^ oo), (21) 

which is equal to the maximally mixed state of two qubits and, of course, not entangled. 

Once more, we see that the same infinite sequence of data leads to radically different 
conclusions for the two priors. 

CONCLUSION 

The most general way in quantum mechanics for obtaining a quantum state from data 
is via the Kraus rule Q. It is clear that the quantum state obtained in this way gener- 
ally depends on some prior state in addition to the data. In an earlier paper |[9|] we have 
established the general principle that a quantum state is never determined by measure- 
ment data alone. This is true even in state preparation, because the prepared state always 
depends on the prior quantum state of the preparation device [9] . 

What we have illustrated here it that this general principle continues to hold in sit- 
uations where measurements are repeated many times. A quantum state is never deter- 
mined by measurement data alone, even in the limit of infinitely many measurements. 
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